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Abstract — One of the unexpected breakdowns in the existing 
theory of quantum serial turbo coding is that a quantum 
convolutional encoder cannot simultaneously be recursive and 
non-catastrophic. These properties are essential for a quantum 
turbo code to have an unbounded minimum distance and for its 
iterative decoding algorithm to converge, respectively. Here, we 
show that the entanglement-assisted paradigm gives a theoretical 
and practical "turbo boost" to these codes, in the sense that 
an entanglement-assisted quantum (EAQ) convolutional encoder 
can possess both of the aforementioned desirable properties, and 
simulation results indicate that entanglement-assisted turbo codes 
can operate reliably in a noise regime 5.5 dB beyond that of 
standard quantum turbo codes. Entanglement is the resource that 
enables a convolutional encoder to satisfy both properties because 
an encoder acting on only information qubits, classical bits, 
gauge qubits, and ancilla qubits cannot simultaneously satisfy 
them. We give several examples of EAQ convolutional encoders 
that are both recursive and non-catastrophic and detail their 
relevant parameters. Finally, simulation results demonstrate that 
interleaved serial concatenation of EAQ convolutional encoders 
leads to a powerful code construction with excellent performance 
on a memoryless depolarizing channel. 

Index Terms — quantum communication, entanglement-assisted 
quantum turbo code, entanglement-assisted quantum error cor- 
rection, recursive, non-catastrophic, entanglement-assisted quan- 
tum convolutional code 



I. Introduction 

Classical turbo codes represent one of the great successes of 
the modern coding era 0, 0, 0, 0. These near Shannon- 
limit codes have efficient encodings, they offer astounding 
performance on memoryless channels, and their iterative de- 
coding algorithm quickly converges to an accurate error esti- 
mate. They are probabilistic codes, meaning that they possess 
sufficient structure to ensure efficient encoding and decoding, 
yet they have enough randomness to allow for analysis of their 
performance with the probabilistic method 0, 0, 0. 

The theory of quantum turbo codes is much younger than its 
classical counterpart 0, 0, and we still stand to learn more 
regarding these codes' performance and structure. Poulin et 
al. set this theory on a firm foundation [7] in an attempt to 
construct explicit quantum codes that come close to achieving 
the quantum capacity of a quantum channel 0,0, ifTOl . fiTTl . 
The structure of a quantum serial turbo code is similar to 
its classical counterpart — one quantum convolutional encoder 
(El, lU3l followed by a quantum interleaver and another 
quantum convolutional encoder. The encoder "closer to the 
channel" is the inner encoder, and the one "farther from the 
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channel" is the outer encoder. One of the insights of Poulin et 
al. in Ref. was to "quantize" the classical notion of a state 
diagram iTBll . 031 — this diagram helps in analyzing important 
properties of the constituent quantum convolutional encoders 
that directly affect the performance of the resulting quantum 
turbo code. 

Despite Poulin et al.'s success in providing a solid the- 
oretical construction, the engines of their quantum turbo 
codes had an unexpected breakdown^] They found that quan- 
tum convolutional encoders cannot be simultaneously non- 
catastrophic and recursive, two desirable properties that can 
hold simultaneously for classical convolutional encoders and 
are the reason underpinning the high performance of classical 
turbo codes 0, 0. These two respective properties ensure 
that an iterative decoder performs well in estimating errors 
and that the turbo code has an unbounded minimum distance 
growing almost linearly with the length of the code 0, 
fT6lL ifTTl . Quantum convolutional encoders cannot have these 
properties simultaneously, essentially because stabilizer op- 
erators must satisfy stringent commutativity constraints in 
order to form a valid quantum code (see Theorem 1 of 
Ref. 0). Thus, the existing quantum turbo codes with non- 
catastrophic constituent quantum convolutional encoders do 
not have unbounded minimum distance, but Poulin et al. con- 
ducted numerical simulations and showed that performance of 
their quantum turbo codes appears to be good in practice. 

The breakdown in the quantum turbo coding theory has led 
researchers to ponder if some modification of the quantum 
turbo code construction could have unbounded minimum 
distance with the iterative decoding algorithm converging fT8l . 
One possibility is simply to change the paradigm for quantum 
error correction, by allowing the sender and the receiver access 
to shared entanglement before communication begins. This 
paradigm is known as the "entanglement-assisted" setting, 
and it simplifies both the theory of quantum error correction 
E2), l20l and the theory of quantum channels |2T1L l22l . In 
entanglement-assisted quantum (EAQ) error correction, it is 
not necessary for a set of stabilizer operators to satisfy the 
stringent commutativity constraints that a standard quantum 
code should satisfy, allowing us to produce EAQ codes from 
arbitrary classical codes |^] In quantum channel theory, the 
classical and quantum capacities of a quantum channel assisted 
by entanglement are the only known capacities for which we 
can claim a complete understanding in the general case — 
both are formulas involving an optimization of the quantum 
mutual information with respect to all pure, bipartite entangled 
inputs to a single use of the channel (formally analogous to 

1 Please excuse the unavoidable pun! 

2 This holds for EAQ convolutional codes in addition to EAQ block 
codes (23], El. 
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Shannon's formula for the classical capacity of a classical 
channel l25ll ). The simplification that entanglement offers to 
the theory of quantum information has led to the following 
musing of Hayden et al. 



"To what extent does the addition of free entangle- 
ment make quantum information theory similar to 
classical information theory?" 

In this paper, we show that entanglement can "boost the 
engines of quantum turbo codes." Specifically, we find many 
examples of EAQ convolutional encoders that can simul- 
taneously be recursive and non-catastrophic. A "quantized" 
version of the result in Ref. |5 ] then implies that a quantum 
serial turbo code employing such an encoder along with 
another non-catastrophic encoder has an unbounded minimum 
distance fTTlL and non-catastrophicity implies that it has good 
iterative decoding performance. This result for EAQ convolu- 
tional encoders holds partly because all four Pauli operators 
acting on half of a Bell state are distinguishable when per- 
forming a measurement on both qubits in the Bell state (much 
like the super-dense coding effect (23). The ability of EAQ 
convolutional encoders to be simultaneously non-catastrophic 
and recursive is another way in which shared entanglement 
aids in a straightforward quantization of a classical result — this 
assumption thus simplifies and enhances the existing theory of 
quantum turbo coding. 

Our simulation results indicate that particular constructions 
of entanglement-assisted turbo codes have excellent perfor- 
mance on a memory less depolarizing channel. In particular, 
we found two quantum convolutional encoders that are com- 
parable to the first and third encoders of Poulin et al. Q, 
in the sense that they have the same number of memory 
qubits, information qubits, ancilla qubits, and a comparable 
distance spectrum. These encoders are non-catastrophic, and 
they become recursive after replacing all of the ancilla qubits 
with ebits. Additionally, the encoders with full entanglement 
assistance have a distance spectrum much improved over 
the unassisted encoders, essentially because entanglement 
increases the ability of a code to correct errors l28l . We 
constructed a quantum serial turbo code with these encoders 
and conducted four types of simulations: the first with the 
unassisted encoders, a second with full entanglement assis- 
tance, a third with the inner encoder assisted, and a fourth 
with the outer encoder assisted. The unassisted quantum turbo 
codes have performance similar to those of Poulin et al. in 
Ref. [7]. The encoders with full entanglement assistance 
have a tremendous boost in performance over the unassisted 
ones, in the sense that they can operate reliably in a noise 
regime several dB beyond the unassisted turbo codes. This 
boost is due to the improvement in the distance spectrum 
and is also due to the encoder becoming recursive. Also, 
these codes come close to achieving the entanglement-assisted 
hashing bound I2T1 . 1291 , which is the ultimate limit on their 
performance. The quantum turbo codes with inner encoder 
entanglement assistance have performance a few dB below 
the fully-assisted code, but one advantage of them is that 
other simulations indicate that they are more tolerant to noise 
on the ebits. The quantum turbo codes with outer encoder 



entanglement assistance have performance not far beyond 
that of the unassisted turbo codes. The hope for this last 
construction was that assisted outer encoders with an improved 
distance spectrum combined with unassisted inner encoders 
satisfying the weaker condition of quasi-recursiveness would 
be sufficient to have a marked improvement over the unassisted 
turbo codes, but our simulation results have shown that this 
intuition does not hold. 

We organize this paper as follows. The first section estab- 
lishes notations and definitions for EAQ codes similar to those 
in Ref. [7 ], and it also shows a way in which an EAQ code with 
only ebits is remarkably similar to a classical code. Section III 
defines the state diagram of an EAQ convolutional encoder — 
it reveals if an encoder is non-catastrophic and recursive, and 



we review how to check for these properties. Section IV gives 
several examples of non-catastrophic, recursive EAQ convolu- 
tional encoders and details their distance spectra. We discuss 
the construction of an EAQ serial turbo code in Section [V] 
and give several combinations of serial concatenations that 
have good minimum-distance scaling. Section [Vl] contains 



our simulation results with accompanying interpretations of 



them. In Section VII we show that entanglement is in fact 
the resource that enables a convolutional encoder to be both 
recursive and non-catastrophic — a corollary of Theorem 1 in 
Ref. Q states that other resources such as classical bits, gauge 



qubits, and ancilla qubits do not help. Section |VIII| then dis- 
cusses encoders that act on information qubits, ancilla qubits, 
ebits, and classical bits, and it gives an example of an encoders 
that can be recursive and non-catastrophic. Section |IX| states 



some general observations regarding recursiveness and non- 
catastrophicity for different types of encoders. The conclusion 
summarizes our contribution and states many open questions. 

II. EAQ Codes 

We first review some important ideas from the theory of 
EAQ codes in order to prepare us for defining convolutional 
versions of them along with their corresponding state dia- 
grams. The development is similar to that in Section III of 
Ref. Q. 

The encoder V of an EAQ code produces an encoded state 
by acting on a set of k information qubits in a state a 
ancilla qubits, and c ebits: 



AB 



\*t) AB ), 



where 



|O a ) A = |O) 0a , 



\*+) AB = j=(\00) AB + \ll) AB ) 



V2 

n = k + a + c, 

and the sender Alice possesses c halves of the entangled 
pairs while the receiver Bob possesses the other c halves 
(see Figure [T]). In what follows, we abuse notation by having 
V refer to the "Clifford group" unitary operator that acts as 
above, but having it also refer to a binary matrix that acts on 
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Fig. 1 . The encoder of an EAQ code. Alice acts on her logical qubits L , local 
ancillas S, and her halves of the ebits E with an encoding unitary V . This 
encoding produces physical qubits P that she then inputs to a noisy quantum 
channel. The entanglement-assisted paradigm assumes that noise does not 
occur on Bob's half of the ebits, but we later study this setting in some 
simulations because it could occur in practice. 



binary vectors — these binary vectors represent different Pauli 
operators that are part of the specification of an EAQ code 
(e.g., see Ref. l20l^ 

Suppose now that Alice transmits her n qubits of the 
encoded state 1^) over a noisy Pauli channel. Then the 

resulting state is P A where P A is some n-fold tensor 

product of Pauli operators (in what follows, we simply say 
an "n-qubit Pauli operator"). Suppose that Bob applies the 



inverse (V A )^ of the encoding to the state P A 
resulting state has the following form: 

(v A yp A \^) AB = (y A yp A v A m A ® \o a ) A 
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L A \iP) A ®S A \O a ) A ®E A \<S>t) AB ), 



where L A is some /c-qubit Pauli operator, S A 
is an a-qubit Pauli operator, and E A = E A (g) • 
c-qubit Pauli operator. Observe that 



(1) 



is some 



s A K) A = \siY 



®\SaY 



where Si = if Si G {/, Z} and Si = 1 otherwise. The fact 
that Si is invariant under the application of a Pauli Z operator 
implies that a quantum code can be degenerate (where two 
different errors mapping to the same syndrome have the same 
effect on the encoded state). Also, observe that 



E A \$t) AB = \$(ei, x ,e 1 , z )) 



AB 



\®(e c ,x,ec,z)) 



AB 



where |$ (e^ x ,e iy 
states and 



AB 

)) denotes the four distinguishable Bell 



(ei, x ,e i>z ) = (0,0) if E A = I, 

(e it x,ei, z ) = (0,l) XE A = Z, 

(e i}X ,e itZ ) = (l,0) ifE A = X, 

(ei, x ,e i>z ) = (l,l) if E A = Y. 

3 The representation of the encoder as a binary matrix leads to a loss 
of global phase information. Though, this global phase information is not 
important because measurement of the syndrome destroys it, and it is not 
necessary for faithful recovery of the encoded state. 



Thus all four Pauli operators are distinguishable in this case 
by performing a Bell measurement (similar to the super-dense 
coding effect l27l ). Ebits do not contribute to the degeneracy 
of a quantum code because different errors E A lead to distinct 
measurement results. 

Bob can perform Z basis measurements on the ancillas and 
Bell measurements on the ebits to determine the syndrome 
r (P) of the Pauli error P: 

Consider the following relation between the binary represen- 
tations of P, V, L, S, and E in {j}: 



PV~ 



(L:S:E). 



The syndrome r (P) only partially determines S, but it fully 
determines E. Let us decompose the binary representation of 
S as S = S z + S x and that of E as E = E z + E x . When 
Bob performs his measurements, he determines S x , E z , and 
E x . That is, he determines the following relations between the 
components Sf of S x and the components of the syndrome 
r(P): 

Sf = X if si = 1, Sf = I otherwise. 

The syndrome also determines the i th components Ef and Ef 
ofE: 



E x 
Ef 



X if e { , 
Z if e ? 2 



= 1, 
1, 



Ef 
Ef 



- I otherwise, 
I otherwise. 



The phenomenon of degeneracy represents the most radical 
departure of quantum coding from classical coding 1301 , 13B . 
Consider two different physical errors P and P' that differ 
only by Z operators acting on the ancillas: 



P 
P' 



(L : S x +S Z : E x + E Z )V 
(L:S X +S /Z :E X + E Z )V 
P + (4 : S z + S' z : I C )V, 



where Ik is a length 2k zero vector (the binary representation 
of a /c-fold tensor product of identity operators). These differ- 
ent errors lead to the same error syndrome. In the classical 
world, this would present a problem for error correction. 
But this situation does not cause a problem in the quantum 
world for the errors P and P' — the logical error affecting the 
encoded quantum information is the same for both P and P' , 
and Bob can correct either of these errors simply by applying 
L~ x after decoding. 

We now define several sets of operators that are important 
for determining the properties and performance of EAQ codes. 
The set C (J) of harmless, undetected errors consists of all 
operators P that have a zero syndrome, yet have no effect on 
the encoded state: 

C(I) = {P:P= (J fc : S z : I C )V, S z G {/, Z} a } . 

This set of operators is equivalent to the isotropic subgroup, in 
the language of Refs. fT9lL (20). It is also analogous to the all- 
zero codeword of a classical code. The set C (L) of harmful, 
undetected errors consists of all operators P that have a zero 
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syndrome, yet change the quantum information in the encoded 
state: 

C(L) = {P:P=(L:S Z :I C )V, S z e{I,Z} a }, 

where L ^ Ik. This set of operators corresponds to a single 
logical transformation on the encoded state, depending on the 
choice of L, and it is analogous to a single codeword of a 
classical code. The set C (L,S X ,E X ,E Z ) corresponds to a 
particular logical transformation and syndrome, and is thus a 
logical coset: 

C{L,S X ,E X ,E Z ) 
= {P : P = (L : S z + S x : E z + E X )V, S z e {/, Z} a } . 

It is analogous to a single erred codeword of a classical code 
if S x , E x , or is non-zero. The operator codewords of an 
EAQ code belong to the following set C: 

C = \JC(L), (2) 

L 

where L is an arbitrary /c-qubit Pauli operator. The set C is 
equivalent to the full set of logical operators for the code, 
and it is analogous to the set of all codewords of a classical 
code. These definitions lead to the definition of the minimum 
distance of an EAQ code as the minimum weight of an 
operator P in C — C (I): 

d (C) = min {w (P) : P e C - C (I)} . 

This definition is similar to the definition of the distance of a 
classical code, but it incorporates the coset structure of a quan- 
tum code. Thus, we can determine the performance of the code 
in terms of distance by tracking its logical operators, and this 
intuition is important when we move on to EAQ convolutional 
codes. Also, the logical operators play an important role in 
decoding because a maximum likelihood decoding algorithm 
for a quantum code estimates the most likely logical error 
given the syndrome and a particular physical noise model^] 

We end this section with a final remark concerning EAQ 
codes and the musing of Hay den et al. in Ref. (26). Suppose 
that an EAQ code does not exploit any ancilla qubits and uses 
only ebits. Then the features of the code become remarkably 
similar to that of a classical code. Degeneracy, a uniquely 
quantum feature of a code, does not occur in this case because 
the syndrome completely determines the error, in analogy 
with error correction in the classical world. Also, the code 
loses its coset structure, so that C (I) is equal to the identity 
operator, C (L) is equivalent to just one logical operator, and 
the definition of the code's minimum distance is the same as 
the classical definition. 

III. EAQ CONVOLUTIONAL CODES 

An EAQ convolutional code is a particular type of EAQ 
code that has a convolutional structure. In past work on 
this topic (231, EH, E3, we adopted the "Grassl-Rotteler" 
approach to this theory (33), by beginning with a mathematical 

4 The definitions for the maximum likelihood decoder of an EAQ code are 
nearly identical to those for stabilizer codes in Section IIIC of Ref. Q. Thus, 
we do not give them here. 



description of the code and determining a Grassl-Rotteler 
pearl-necklace encoder that can encode it. Here, we adopt 
the approach of Poulin et al. (7) which in turn heavily 
borrows from ideas in classical convolutional coding (14). We 
begin with a seed transformation (a "convolutional encoder" 
or a "quantum shift-register circuit") and determine its state 
diagram, which yields important properties of the encoder. We 
can always rearrange a Grassl-Rotteler pearl-necklace encoder 
as a convolutional encoder (34), (35) , (36], but it is not clear 
that every convolutional encoder admits a form as a Grassl- 
Rotteler pearl necklace encoder. For this reason and others, 
we adopt the Poulin et al. approach in what follows. 

An EAQ convolutional encoder is a "Clifford group" unitary 
U that acts on m memory qubits, k information qubits, a 
ancilla qubits, and c halves of ebits to produce a set of m 
memory qubits and n physical or channel qubits J^] where 
n = k + a + c. The transformation that it induces on binary 
representations of Pauli operators acting on these registers is 
as follows: 

(M' : P) = (M :L:S : E)U, 

where M' acts on the m output memory qubits, P acts on 
the n output physical qubits, M acts on the m input memory 
qubits, L acts on the k information qubits, S acts on the a 
ancilla qubits, and E acts on the c halves of ebits. Although 
the quantum states in these registers can be continuous in 
nature, the act of syndrome measurement discretizes the errors 
acting on them, and the above classical representation is useful 
for analysis of the code's properties and the flow of the 
logical operators through the encoding circuit (recall that the 
goal of a decoding algorithm is to produce good estimates 
of logical errors). This representation is similar to the shift- 
register representation of classical convolutional codes, with 
the difference that the representation there corresponds to the 
actual flow of bit information through a convolutional circuit, 
while here it is merely a useful tool for analyzing the properties 
of the encoder. 

The overall encoding operation for the code is the trans- 
formation induced by repeated application of the above seed 
transformation to a quantum data stream broken up into 
periodic blocks of information qubits, ancilla qubits, and 
halves of ebits while feeding the output memory qubits of 
one transformation as the input memory qubits of the next 
(see Figure 6 of Ref. [7 ] for a visual aid). The advantage of a 
quantum convolutional encoder is that the complexity of the 
overall encoding scales only linearly with the length of the 
code for a fixed memory size, while the decoding complexity 
scales linearly with the length of the code by employing a local 
maximum likelihood decoder combined with a belief propa- 
gation algorithm [7 ]. The quantum communication rate of the 
code is essentially k/n while the entanglement consumption 
rate is c/n, if the length of the code becomes large compared 
to n. 

5 Physical or channel qubits in the entanglement-assisted paradigm are the 
ones that Alice transmits over the channel. 



5 




Fig. 2. (Color online) Seed transformation for an n = 2, k = 1, c = 1, 
m = 1 EAQ convolutional code. 



A. State Diagram 

The state diagram of an EAQ convolutional encoder is the 
most important tool for analyzing its properties, and it is 
the formal quantization of a classical convolutional encoder's 
state diagram 031, El - It examines the flow of the logical 
operators through the encoder with a finite- state machine 
approach, and this representation is important for analyzing 
both its distance and its performance under the iterative 
decoding algorithm of Ref. Q. The state diagram is a directed 
multigraph with 4 m vertices that we think of as memory 
states. We label each memory state with an m-qubit Pauli 
operator M. We connect two vertices with a directed edge 
from M — ^ M' , labeled as (L,P), if there exists a /c-qubit 
Pauli operator L, an n-qubit Pauli operator P, and an a-qubit 
Pauli operator S z e {I, Z} a such that 

(M' : P) = (M : L : S z : J C )E7. (3) 

We refer to the labels L and P of an edge as the respective 
logical and physical label. 

As an example, consider the transformation depicted in 
Figure [2] It acts on one memory qubit, one information qubit, 
and one half of an ebit to produce two channel qubits and 
one memory qubit. Figure [3] illustrates the state diagram 
corresponding to this transformation. There are four memory 
states because there is only one memory qubit, and there are 
16 edges because there are four memory states and four logical 
operators for one information qubit and one ebit. 

B. Non-catastrophicity 

We now recall the definition of non-catastrophicity in 
Ref. 0, which is the formal quantization of the definition 
in Section IX of Ref. lfT4l . The definition from Ref. (7) is the 
same for an EAQ convolutional encoder because it depends 
on the iterative decoding algorithm used to decode the code, 
and we can exploit the same iterative decoding algorithm as in 
Ref. Q to decode EAQ convolutional codes. A path through 
the state diagram is a sequence Mi , . . . , M t of vertices such 
that Mi — )> Mi+i is an edge belonging to it. Each logical 
operator of C in ^ corresponds to a path in the state diagram, 
with the sequence of vertices in the path being the states of 
memory traversed while encoding the logical operator. The 
physical and logical weights of a logical operator are equal to 
the sums of the corresponding weights of the edges traversed 
in a path that encodes the logical operator. A zero physical- 
weight cycle is a cycle in the state diagram such that all 
edges in the cycle have zero physical weight. Finally, an EAQ 
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Fig. 3. The state diagram corresponding to the seed transformation in 
Figure [2] The state diagram allows us to check that the encoder is non- 
catastrophic and non-recursive. 



encoder acting on m memory qubits, k information qubits, a 
ancilla qubits, and c ebits is non-catastrophic if every zero 
physical-weight cycle in its state diagram has zero logical 
weight. 

Why is this definition an appropriate definition of non- 
catastrophicity? We briefly review the argument in Ref. 171 . 
First, consider Figure 8 of Ref. |7|. There, the seed transfor- 
mation is the same as that in our Figure [2j with the exception 
that the ebit is replaced by an ancilla qubit. Such a replacement 
leads to a doubling of the number of edges in the state diagram 
in our Figure [3] because the state diagram should include all of 
the transitions where a Z operator acts on the ancilla qubits 
(these all lead to other logical operators in their case). The 
new state diagram (see Figure 9 of Ref. |7]) includes a self- 
loop at memory state Z with zero physical weight and non- 
zero logical weight, and it is thus catastrophic according to 
the definition. The problem with such a loop is that it can 
"throw off the iterative decoding algorithm. Consider again 
Figure 8 of Ref. Q. If a Y error were to act on the second 
physical qubit in one frame of the code (while the identity 
acts on the rest), then "pushing" this error through the inverse 
of the encoder applies an X operator to the ancilla qubit and 
produces one syndrome bit for that frame of the code, but 
it propagates Z errors onto every logical qubit in the stream 
while applying Z errors to every memory qubit. All of these 
other errors go undetected by an iterative decoder because the 
error propagation does not trigger any additional syndrome 
bits besides the initial one that the Y error triggered. 

Observe that the state diagram in Figure [3] for the EAQ 
convolutional encoder does not feature a zero physical-weight 
cycle with non-zero logical weight. Thus, the encoder is 
non-catastrophic, illustrating another departure from classical 
theory. Non-catastrophicity in the quantum world is not only 
a property of the encoder, but it also depends on the resources 
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available for encoding. If we analyze the above scenario 
that leads to catastrophic error propagation in the encoder of 
Figure 8 of Ref. Q, we see that it does not lead to such 
propagation for our setting in Figure [2] Indeed, suppose again 
that a Y error acts on the second physical qubit in a particular 
frame of the code (with the identity acting on all other qubits). 
Pushing this Y error through the inverse of the encoder leads to 
an X operator acting on the ebit and a Z operator acting on the 
memory. The X operator is detectable by a Bell measurement, 
and the Z operator acting on the memory propagates to a 
Z operator acting on an ebit in the next frame and then to 
all ebits and information qubits in successive frames. These 
Z operators acting on the ebits are all detectable by a Bell 
measurement (whereas they are not detecable when acting on 
an ancilla), so that the iterative decoding algorithm can still 
correct for errors in the propagation because all these errors 
trigger syndromes. 

C. Recursiveness 

Recursiveness is another desirable property for an EAQ con- 
volutional encoder when it is employed as the inner encoder 
of a quantum serial turbo code|^Tfiis property ensures that a 
quantum serial turbo code on average has a distance growing 
near-linearly with the length of the code, and the proof of this 
statement (T6), ifTTl is essentially a direct translation of the 
classical proof |5|. In short, an EAQ convolutional encoder 
is quasi-recursive if it transforms every weight-one Pauli 
operator JQ, 1^, or Zi to an infinite- weight Pauli operator Q. 
It is recursive if, in addition to being quasi-recursive, every 
element in the logical cosets C PQ), C (Yj), and C (Zi) 
has infinite weight. This stringent requirement is necessary 
due to the coset structure of EAQ codes Q One might think 
that quasi-recursiveness is a sufficient requirement for good 



performance, but our simulation results in Section [VI] indicate 
that recursiveness is indeed necessary because a turbo code 
has a significant boost in performance if its inner encoder is 
recursive. 

It seems like it would be demanding to determine if this 
condition holds for every possible input, but we can exploit 
the state diagram to check it. The algorithm for checking 
recursiveness is as follows. First, we define an admis sable path 
to be a path in which its first edge is not part of a zero physical- 
weight cycle Q. Now, consider any vertex belonging to a 
zero-physical weight loop and any admissable path beginning 
at this vertex with logical weight one. The encoder is recursive 
if all such paths do not contain a zero physical- weight loop Q. 
The idea is that a weight-one logical operator is not enough 
to drive a recursive encoder to a memory state that is part of a 
zero-physical weight loop — the minimum weight of a logical 
operator that does so is two. 

The example encoder in Figure [2] is not recursive. The only 
vertex in the state diagram in Figure [3] belonging to a zero 
physical-weight cycle is the vertex labeled /. If the logical 

6 Recall that the inner encoder is the one closer to the channel, and the 
outer encoder is the one farther from the channel. 

7 This definition thus implies that a recursive encoder is infinite-depth (it 
transforms a finite- weight Pauli operator to an infinite-weight one) [33], (23], 
but the other implication does not necessarily have to hold. 



input to the encoder is a Z operator followed by an infinite 
sequence of identity operators, the circuit outputs ZZ as the 
physical output, return to the memory state /, and then output 
the identity for the rest of time. Thus, the response to this 
weight-one input is finite, and the circuit is non-recursive. 



Although this example is non-recursive, Section IV details 
many examples of EAQ convolutional encoders that are both 
recursive and non-catastrophic. 

D. Distance Spectrum 

We end this section by reviewing the performance measures 
from Ref. Q, which are also quantizations of the classical 
measures [14]. The distance spectrum F (w) of an EAQ 
convolutional encoder is the number of admissable paths 
beginning and ending in memory states that are part of a 
zero physical-weight cycle, where the physical weight of each 
admissable path is w and the logical weight is greater than 
zero. This performance measure is most similar to the weight 
enumerator polynomial of a quantum block code (37), which 
helps give an upper bound on the error probability of a non- 
degenerate quantum code on a depolarizing channel under 
maximum likelihood decoding [38]. The distance spectrum 
incorporates the translational invariance of a quantum con- 
volutional code and gives an indication of its performance on 
a memoryless depolarizing channel. The next section details 
a simple method to compute the distance spectrum by using 
the state diagram and ideas rooted in Refs. lfT4l . l39lL l40lL 
flTTl . The free distance of an EAQ convolutional encoder is the 
smallest weight w for which F (w) > 0, and this parameter 
is one indicator of the performance of a quantum serial turbo 
code employing constituent convolutional encoders. Although 
one of the applications of our EAQ convolutional encoders 
are as the inner encoders in a quantum serial turbo coding 
scheme, we can also have them as outer encoders in a quantum 
serial turbo code and use the free distance to show that its 
minimum distance grows near-linearly when combined with a 
non-catastrophic, recursive inner encoder. 

IV. Example Encoders 

Our first example of an EAQ convolutional encoder is the 
simplest example that is both recursive and non-catastrophic. 
It exploits one memory qubit and one ebit to encode one 
information qubit per frame. We discuss this first example in 
detail, verifying its non-catastrophicity and recursiveness, and 
we give a method to compute its distance spectrum. We then 
give tables that detail our other examples. We found all of these 
examples by picking encoders uniformly at random from the 
Clifford group, according to the algorithm in Section VI-A.2 
of Ref. 1421 . 

The seed transformation for our first example is as follows: 
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I 
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X 

I 
I 



I 
z 
I 
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X 

I 



I 
I 
z 
I 
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X 
Y 
Z 



Z I 
X z 
X Y 
XXX' 

Y I Y 

Y X Y 



(4) 
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where the first input qubit is the memory qubit, the second 
input qubit is the information qubit, the third is Alice's half 
of the ebit, the first output qubit is the memory qubit, and the 
last two outputs are the physical qubits. We can abbreviate 
the above encoding by taking the binary representation of 
each row at the output and encoding it as a decimal number. 
For example, the first row ZIX has the binary representation 
[ 1 | 1 ] which is the decimal number 33. 
Thus, we can specify this encoder as {33, 29, 30, 7, 45, 47}(^ 
The seed transformation in ^ leads to the state diagram 
of Figure |4| by exploiting ([3]). We can readily check that 
the encoder is non-catastrophic and recursive by inspecting 
Figure [4] The only cycle with zero physical weight is the self- 
loop at the identity memory state with zero logical weight. The 
encoder is thus non-catastrophic. To verify recursiveness, note 
again that the only vertex belonging to a zero physical-weight 
cycle is the self-loop at the identity memory state. We now 
consider all weight-one admissable paths that begin in this 
state. If we input a logical X, we follow the edge (X, IY) 
to the Y memory state. Inputting the identity operator for the 
rest of time keeps us in the self-loop at memory state Y while 
still outputting non-zero physical weight operators. We can 
then check that inputting a Z logical operator followed by 
identities keeps us in the self-loop at the X memory state 
while still outputting non-zero physical weight operators, and 
inputting a Y logical operator followed by identities keeps us 
in the self-loop at the Z memory state. The encoder is thus 
recursive. 

There is a straightforward way to compute the distance 
spectrum of the encoder. This technique borrows from similar 
ideas in the classical theory of convolutional coding lfT4l . 
l39lL (401, ED- We would like to know the number of 
admissible paths with a particular weight beginning and ending 
in memory states that are part of a zero physical- weight cycle. 
Again, for our example, the identity memory state is the only 
memory state part of a zero physical-weight cycle. We create 
a weight adjacency matrix whose entries correspond to edges 
in the state diagram. This matrix has x w in entry (i, j) if there 
is a phsyical-weight-w edge from vertex i to vertex j (with 
the exception of the self-loop at the identity memory state). 
The weight adjacency matrix A for our example is 



A = 



'0 


x 2 


X 


X 


x 2 


x 2 


x 2 


x 2 


x z 


X 


X 


2 

x A 


x 2 


X 


x 1 


X 



where the ordering of vertices is /, X, Y, and Z. Note that 
we place a zero in the (1,1) entry because we do not want to 
overcount the number of admissable paths starting and ending 
in memory states that are part of a zero physical- weight cycle. 
If we would like the number of admissable paths up to an 
arbitrary weight that start and end in the identity memory state, 
then we compute the (1,1) entry of the following matrix: 

(I - A)' 1 = / + A + A 2 + A 3 + -- -. 

The coefficient of x w in the polynomial entry B (1, 1) is the 
number of admissable paths with weight w starting and ending 

8 Note that this convention is different from that of Poulin et al. in Ref. Q- 




(/JO) aw 



Fig. 4. The state diagram corresponding to the seed transformation in (|4j, 
acting on one memory qubit, one information qubit, one ebit, and producing 
two physical qubits. The above state diagram allows us to check that the 
encoder is both recursive and non- catastrophic. 

in the identity memory state. We can approximate the distance 
spectrum, e.g., by computing the matrix B where 

10 

B = $>\ 

i=i 

Computing the above matrix can be computationally expen- 
sive, but we can dramatically reduce the number of computa- 
tions by truncating the polynomial entries of A above degree 
ten before performing each multiplication. For our example, 
the first ten entries of the distance spectrum polynomial are 

2x 3 + 5x 4 + 6x 5 + 23x 6 + 54x 7 + I22x 8 + 29Sx 9 + 737x 10 . 

These coefficients appear in the second column of Table [TT] 
as the first ten values of the distance spectrum for this 
first example encoder. Note that there are faster ways of 
computing the distance spectrum for classical convolutional 
codes [43], and it remains open to determine how to exploit 
these techniques for quantum convolutional encoders. 

Table |T| lists the specifications of many other examples of 
EAQ encoders that are both recursive and non-catastrophic — 
a computer program helped check that these properties hold 
for each of the examples l44l . Included in the list of example 
encoders are some which act on ancilla qubits in addition to 
ebits. These examples demonstrate that we do not necessarily 
require the auxiliary resource of an EAQ convolutional en- 
coder to be ebits alone in order for the encoder to possess 
both properties. Table [TT] gives the distance spectrum for each 
of these encoders. 

V. EAQ Turbo Codes 

We comment briefly on the construction and minimum dis- 
tance of a quantum serial turbo code that employs our example 
encoders as a constituent encoder (the next section details our 
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TABLE I 

Specification of our example EAQ convolutional encoders that are recursive and non-catastrophic. The first column indexes 

THE DIFFERENT ENCODERS, AND ALL OF THE OTHER COLUMNS CORRESPOND TO A PARTICULAR ENCODER. COLUMN "M" GIVES THE NUMBER OF 
MEMORY QUBITS, COLUMN "L" GIVES THE NUMBER OF INFORMATION QUBITS, COLUMN "A" GIVES THE NUMBER OF ANCILLA QUBITS, AND COLUMN 
"E" GIVES THE NUMBER OF EBITS. THE COLUMN LABELED "SEED TRANSFORMATION" SPECIFIES THE SEED TRANSFORMATION FOR EACH ENCODER, 

USING A DECIMAL REPRESENTATION. THE CONVENTION FOR QUBIT ORDERING IS THAT AN ENCODER ACTS ON MEMORY QUBITS, INFORMATION 
QUBITS, ANCILLA QUBITS, AND EBITS TO PRODUCE OUTPUT MEMORY QUBITS AND PHYSICAL OR CHANNEL QUBITS. THE CONVENTION FOR THE BIT 
REPRESENTATION OF A TENSOR PRODUCT OF PAULI OPERATORS IS ZXYZ MAPS TO [1011 1 01 10], AND THE DECIMAL REPRESENTATION OF THIS BIT 
VECTOR IS 182 (SEE THE DISCUSSION AROUND ([4) IN THE TEXT). THE LAST COLUMN GIVES THE FREE DISTANCE OF THE ENCODER. 
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TABLE II 

The distance spectrum for the first seven of our example EAQ 

CONVOLUTIONAL ENCODERS. THE FIRST COLUMN ON THE LEFT GIVES 
THE FIRST 1 1 VALUES OF THE WEIGHT W, AND THE OTHER COLUMNS 
GIVE THE FIRST 1 1 VALUES OF THE DISTANCE SPECTRUM F(w) OF THE 
FIRST SEVEN OF OUR EAQ CONVOLUTIONAL ENCODERS FROM TABLe|T] 



simulation results with different encoders). The construction 
of an EAQ serial turbo code is the same as that in Ref. (7) 
(see Figure 10 there), with the exception that we assume that 
Alice and Bob share entanglement in the form of ebits before 
encoding begins. Alice first encodes her stream of information 
qubits with the outer encoder, performs a quantum interleaver 
on all of the qubits, and then encodes the resulting stream with 
the inner encoder. The quantum communication rate of the 
resulting EAQ turbo code is k 0ut /n ln = (k 0ut /n 0ut )(k ln /n ln ) 
where k 0ut is the number of information qubits encoded by 
the outer encoder, n 0ut is the number of physical qubits output 
from the outer encoder, and a similar convention holds for k ln , 
n In , and the inner encoder. In order for the qubits to match 
up properly, n 0ut must be equal to k ln . The entanglement 
consumption rate of the code is (c 0ut + c In )/n In where c 0ut 
and c In are the total number of ebits consumed by the outer 
and inner encoder, respectively. 



Perhaps the best combination for an EAQ serial turbo code 
is to choose the inner quantum convolutional encoder to 
be a recursive, non-catastrophic EAQ convolutional encoder 
and the outer quantum convolutional encoder to be a non- 
recursive, non-catastrophic standard quantum convolutional 
encoder. This combination reduces entanglement consumption, 
ensures good iterative decoding performance, and ensures that 
the minimum distance of the quantum serial turbo code grows 
as N( d *~ 2 ^ d * where N is the length of the code and d* is the 
free distance of the outer quantum convolutional encoder. The 
proof of this last statement ifTTl . lfT6l is essentially identical 
to the classical proof (5j. Section [VI- A| shows that this com- 
bination also performs well in practice if noise occurs on the 
ebits. Additionally, choosing the outer quantum convolutional 
encoder to encode a highly degenerate code may increase the 
number of errors that the quantum turbo code can correct, in 
a vein similar to the results in Ref. OTI (though we have not 
yet fully investigated this possibility). 

We can construct EAQ serial turbo codes, by serially 
concatenating some of our example "WH encoders" in Table |T| 
with the "PTO encoders" in Table 1 of Ref. |7 ] (ordered from 



left to right). Table [HI] details these different combinations, 
giving their rates and average minimum distance scaling. 



The columns of Table III give the outer encoder, the inner 
encoder, the quantum communication rate, the entanglement 
consumption rate, and the average minimum distance growth 
of a particular EAQ turbo code. The first four combinations 
all have a good average minimum distance scaling, but the 
catalytic rat^] for each of them is zero. The last six combi- 
nations both have a good average minimum distance scaling 
and a positive catalytic rate. 

9 The catalytic rate is the difference between the quantum communication 
rate and the entanglement consumption rate 1 201 . 
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TABLE III 

Various combinations of encoders for quantum turbo codes 
and their minimum distance scaling. 
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VI. Simulation Results 

We performed several simulations of EAQ turbo codes and 
detail the results in this section. This section begins with 
a description of the parameters of the constituent quantum 
convolutional encoders. We then describe how the simulation 
was run, and we finally discuss and interpret the simulation 
results. 

We would like to compare our results with those of Poulin 
et al. Q — our original intent was to exploit their encoders, but 
the lack of a clear exposition of their decimal representation 
convention has excluded us from doing so. So we instead 
randomly generated and filtered encoders that have compa- 
rable attributes to their encoders. The first encoder, dubbed 
"PTOIR," acts on three input memory qubits, one information 
qubit, and two ancillas to produce three output memory qubits 
and three physical qubits. Its decimal representation is 

{1355, 2847, 558, 2107, 3330, 739, 
2009,286,473,1669,1979,189} ' 

and its distance spectrum polynomial is 

llx 5 + 47x 6 + 2b3x 7 + 1187x 8 + 

6024x 9 + 30529x 10 + 153051a 11 + 771650x 12 . 

This encoder is non-catastrophic and quasi-recursive, and its 
parameters and distance spectrum are comparable to those of 
Poulin et al.'s first encoder (though they did not comment 
there whether theirs is quasi-recursive) . Replacing the two 
ancilla qubits in the encoder with two ebits gives a dramatic 
improvement in the distance spectrum polynomial: 

2x 9 + x 10 + bx 11 + Sx 12 + llx 13 + 25x 14 + 

56x 15 + 102x 16 + 217x 17 + 387x 18 + 787x 19 . 

The encoder also becomes recursive after this replacement. Let 
"PTOIREA" denote the EA version of this encoder. 

Our second encoder, dubbed "PT03R," acts on four input 
memory qubits, one information qubit, and one ancilla qubit to 
produce four output memory qubits and two physical qubits. 
Its decimal representation is 

{3683, 3556, 2872, 2211, 3561, 3534, 
729,3136,743,2643,1330,1656} ' 



Fig. 5. The figure plots the hashing bound for quantum codes (solid blue 
curve), the entanglement-assisted hashing bound for EAQ codes (the dash- 
dotted green curve), and the entanglement consumption rate of the father 
protocol that meets the entanglement- assisted hashing bound (dashed red 
curve). The figure plots these different bounds as a function of the depolarizing 
parameter p. Observe that the noise limit for a rate 1/4 quantum code is 
« 0.12689, the noise limit for a rate 1/9 quantum code is « 0.16028, the 
noise limit for a rate 1/4 EAQ code is « 0.35454, and the noise limit for a 
rate 1/9 EAQ code is w 0.49088. We should compare our simulation results 
with these noise limits to observe how close our code constructions come to 
achieving capacity. 



and its distance spectrum polynomial is 

12x 5 + 93x 6 + 600x 7 + 4320x 8 + 31098x 9 + 

224014x 10 + 1604435X 11 + 11469935x 12 . 

The encoder is non-catastrophic and quasi-recursive with its 
parameters and distance spectrum comparable to those of 
Poulin et a/.'s third encoder. Replacing the ancilla with an ebit 
gives an improvement to the distance spectrum polynomial: 

x 6 + 3x 7 + 7x 8 + 29x 9 + 88x 10 + 237 x n + 

716x 12 + 2166x 13 + 6245x 14 + 18696x 15 + 
55889x 16 + 165971x 17 + 492805x 18 + 1465529x 19 , 

and the encoder becomes recursive after this replacement. Let 
"PT03REA" denote the EA version of the "PT03" encoder. 

Our turbo codes consist of interleaved serial concatenation 
of the above encoders with themselves, and we varied the 
auxiliary resources of the encoders to be ebits, ancillas, or 
both. Concatenating PTOIR with itself leads to a rate 1/9 
quantum turbo code, and concatenating PT03R with itself 
leads to a rate 1/4 turbo code. 

We simulated the performance of these EAQ turbo codes 
on a memoryless depolarizing channel with parameter p. The 
action of the channel on a single qubit density operator p is 
as follows: 

p^(l-p) p+ l(XpX + Y P Y + ZpZ). 

The benchmarks for the performance of any standard or 
entanglement-assisted quantum code on a depolarizing channel 
are the hashing bounds (45), (2T). The quantum hashing bound 
determines the rate at which a random quantum code can 
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operate reliably for a particular depolarizing parameter p, but 
note that it does not give the ultimate capacity limit because 
degeneracy can improve performance (3TJ, (46l , (47l . (48l . On 
the other hand, the entanglement- as sis ted hashing bound does 
give the ultimate capacity at which an EAQ code can operate 
reliably for a particular depolarizing parameter p, whenever 
an unbounded amount of entanglement is available (2T1 . (29l . 
The hashing bound for quantum communication is 

l-iH 2 (p)+p\og 2 (3)], 

and the hashing bound for an entanglement-assisted quantum 
code is 

l--[ J ff 2 (p)+plog 2 (3)]. 

The father protocol is a quantum communication protocol 
that achieves the entanglement-assisted quantum capacity (49l , 
1551 , I5T1L while attempting to minimize its entanglement 
consumption rate. Its entanglement consumption rate is prov- 
ably optimal for certain channels such as the dephasing 
channel 1521 , (53l , 1541 , but it is not necessarily optimal for 
the depolarizing channel. The entanglement consumption rate 
of the father protocol on a depolarizing channel is 

1 -[H 2 {p)+p\og 2 (3)]. 

This entanglement consumption rate implies that the only 
resources involved in an encoded transmission of the father 
protocol are information qubits and ebits because the sum 
of its quantum communication rate and its entanglement 
consumption rate is equal to one. Figure [5] plots these different 
hashing bounds and depicts the location of these bounds for 
our different quantum turbo codes. We performed Monte Carlo 
simulations to determine the performance of our example 
EAQ turbo codes when decoded with an iterative decoding 
algorithm (the iterative decoding algorithm is similar to that 
from Ref. |7), with the exception that it decodes EAQ turbo 
codes). We developed a Matlab computer program (44), and 
we conducted these simulations for about one month on 
approximately 70 Ubuntu Linux servers in the School of 
Computer Science at McGill University. 

A single run of a simulation selects a quantum turbo code 
with a particular number of logical qubits and a random choice 
of interleave^ and it then generates a Pauli error randomly 
according to a depolarizing channel with parameter p. This 
Pauli error leads to an error syndrome, and the syndrome 
and channel model act as inputs to the iterative decoding 
algorithm. The iterative decoding algorithm terminates if the 
hard decision on a recovery operation is the same decision 
from the previous iteration, or it terminates after a maximum 
number of iterations (though we never observed the number 
of iterations until convergence exceeding six). One run of 
the simulation declares a failure if the estimated recovery 
operation is different from the correct recovery operation. The 
ratio of simulation failures to the total number of simulation 
runs is the word error rate (WER). We ran every configuration 
(choice of code, depolarizing parameter, and number of logical 
qubits) until we observed at least 100 failures — this number 
gaurantees a reasonable statistical confidence in the results of 
the simulations. 



Our first simulation involved the serial concatenation of 
the unassisted PTOIR encoder with itself, and Figure [6ja) 
displays the results. The performance is similar to that of the 
first encoder of Poulin et al. (see Figure 12 of Ref. 171). A 
noticeable feature of Figure [6ja) is the existence of a pseu- 
dothreshold, such that increasing the number of encoded qubits 
of the turbo code decreases the WER for all depolarizing noise 
rates below the pseudothreshold. The pseudothreshold is not a 
true threshold because this particular turbo code has a bounded 
minimal distance, and at some point, the WER should begin 
increasing if we continue to increase the number of encoded 
qubits [7]. The pseudothreshold occurs at a depolarizing noise 
rate « 0.098, and this pseudothreshold is within 2.14 dB of 
the 0.16028 noise limit for a rate 1/9 code. 

Our second simulation tested the serial concatenation of 
the PTOIREA encoder with itself, and Figure [6jb) displays 
the results. This turbo code uses the maximal amount of 
entanglement at a rate of 8/9 and thus is an instance of the 
father protocol. Entanglement assistance gives the turbo code 
a dramatic boost in performance, in the sense that it can with- 
stand far higher depolarizing noise rates than the unassisted 
turbo code. The true threshold occurs at « 0.345 — it is a true 
threshold because the code has a provably unbounded minimal 
distance ifTTl , and the WER should continue to decrease as we 
increase the number of encoded qubits. This threshold is 5.47 
dB beyond the pseudothreshold of the unassisted turbo code, 
and it is within 1.53 dB of the 0.49088 noise limit from the 
EA hashing bound. Thus, this code construction is operating in 
a noise regime in which standard quantum codes are simply 
not able to operate (compare with the results in Refs. 1551 , 

ESI, (571, 0, EE EH, (551, ED). 

Our third simulation tested the serial concatenation of the 
PTOIREA inner encoder with the PTOIR outer encoder, and 
this EAQ turbo code has an entanglement consumption rate of 
6/9. The inner encoder is recursive, but the outer encoder's 
free distance is not as high as that in the previous simulation. 
Thus, we should expect performance not to be as high as in 
the previous simulation. Figure |6jc) displays the results. Since 
the inner encoder is recursive, the EAQ turbo code has an 
unbounded minimum distance, and the threshold in Figure [6jc) 
is a true threshold. It occurs approximately at a depolarizing 
noise rate of 0.28, which is within 0.91 dB of the previous 
threshold and within 2.44 dB of the 0.49088 noise limit from 
the EA hashing bound. 

Our final simulation tested the serial concatenation of the 
PTOIR inner encoder with the PTOIREA outer encoder, and 
this EAQ turbo code has an entanglement consumption rate 
of 2/9. The inner encoder is quasi-recursive, and the outer 
encoder has a significantly higher free distance than in the 
previous simulation because it has entanglement assistance. 
This simulation provides a good test to determine the effec- 
tiveness of an inner encoder that is quasi-recursive. That is, 
one might think that quasi-recursiveness of the inner encoder 
combined with an outer encoder with high free distance would 
be sufficient to produce a turbo code with good performance 
(the thought is that this would explain the good performance 
in the first simulation), but Figure [6jd) demonstrates that 
this intuition does not hold. The pseudothreshold occurs at 
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Fig. 6. The figure plots the results of four different simulations that exploit the PTOIR and PTOIREA encoders. The resulting EAQ turbo codes have rate 
1/9 and varying entanglement consumption rates. Each plot demonstrates the existence of a threshold, where the code performance increases whenever the 
channel noise rate is below the threshold. The vertical dotted line in each plot indicates the approximate location of the threshold. 



approximately 0.14, but this value is only a rough estimate 
because it is not as clearly defined. This pseudothreshold is 
only 1.55 dB higher than the threshold for the unassisted code, 
and it is 5.45 dB away from noise limit from the EA hashing 
bound. 

We conducted similar simulations with the PT03R and 
PT03REA encoders, and Figures |7ja-d) display the results. 
The entanglement consumption rates in Figures [7Ja-d) are 0, 
3/4, 2/4, and 1/4, respectively. The results are similar to 
the previous simulations with the difference that the noise 
limits and thresholds are lower because these codes have 
higher quantum transmission rates. The thresholds (pseudo 
and true) occur at approximately 0.067, 0.216, 0.17, and 0.09 
in Figures [7Ja-d), respectively. There is one feature of the 
result in Figure [TJc) that is worth discussing. Consider that 
the midpoint of the highest threshold 0.216 and the threshold 



0.09 is 0.153. One might expect then that a turbo code with 
entanglement consumption rate halfway between 1/4 and 3/4 
should have a threshold occurring at the halfway point 0.153. 
But this is not the case for the code in Figure [TJc) with 
entanglement consumption rate 2/4. Its threshold occurs at 
0.17, which is slightly higher than 0.153. Thus, this EAQ turbo 
code is making judicious use of the available entanglement by 
placing the two ebits in the inner encoder. We would have to 
conduct further simulations to determine if placing one ebit 
in the outer encoder and one in the inner encoder would do 
better, but our suspicion is that the aforementioned use of 
entanglement is better because the inner encoder is recursive 
under this choice. 

A. Noise on Ebits 

We conducted another set of simulations to determine how 
noise on Bob's half of the ebits affects the performance of a 
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Fig. 7. The figure plots the results of four different simulations that exploit the PT03R and PTO3REA encoders. The resulting EAQ turbo codes have rate 
1/4 and varying entanglement consumption rates. Each plot demonstrates the existence of a threshold, where the code performance increases whenever the 
channel noise rate is below the threshold. The vertical dotted line in each plot indicates the approximate location of the threshold. 



code. This possibility has long been one of the important prac- 
tical questions concerning the entanglement-assisted paradigm, 
and some researchers have provided partial answers (62), (63) , 
(64), 1 65]. We briefly review some of these contributions. 
Shaw et al. first observed that the [[7,1,3]] Steane code is 
equivalent to a [[6, 1, 3; 1]] entanglement-assisted code that can 
also correct a single error on Bob's half of the ebit. This ob- 
servation goes further: any standard, non-degenerate [[n, /c, d]] 
quantum code is equivalent to an [[n — d + 1, k : d\ d — 1]] 
entanglement-assisted code that can correct any |_(^ — 1) /2J 
errors on Alice and Bob's qubits. This result holds because 
tracing over any n — d + 1 qubits in the original standard 
code gives a maximally mixed state on d — 1 qubits, and 
the purification of these qubits are n — d + 1 encoded halves 
of ebits that Alice possesses 1 66 1. Wilde and Fattal observed 
that entanglement-assisted codes correcting for errors on Bob's 



side slightly improve the threshold for quantum computation 
if ebit errors occur less frequently than gate errors (63) . Hsieh 
and Wilde studied this question in the Shannon-theoretic con- 
text and determined an expression for capacity when channel 
errors and entanglement errors occur [64]. As a side note, Lai 
and Brun recently looked for codes attempting to achieve the 
opposite goal (28) — their codes try to maximize the number 
of channel errors that can be corrected while minimizing the 
correction on Bob's half of the ebits. 

Figure [8] plots the results of our simulations that allow 
for noise on Bob's half of the ebits. We performed three 
different types of simulations: the first was with the PT03REA 
inner / PT03R outer combination, the second with the PT03R 
inner / PT03REA outer combination, and the third with the 
PT03REA inner / PT03REA outer combination. For each 
simulation, we kept the channel noise rates fixed at 0.14, 0.07, 
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Performance with noise on ebits 




Depolarizing rate on ebits 



Fig. 8. The figure plots the performance of different inner and outer encoder 
combinations as a function of noise rate on Bob's half of the ebits. It is 
surprising that the combination PT03REA inner / PT03R outer can withstand 
more ebit noise than the combination PT03R inner / PTO3REA outer even 
though it uses more ebits than this latter combination. 



and 0.19, respectively, because the codes already performed 
reasonably well at these noise rates, and our goal was to 
understand the effect of ebit noise on code performance. The 
codes all performed about the same as they did without ebit 
noise if we set the ebit noise rate at 10 -5 . Increasing the 
ebit noise rate an order of magnitude to 10 -4 has the least 
effect on the PT03REA inner / PT03R outer combination. 
Increasing it further to 10 -3 deteriorates the performance 
of the other combinations while still having the least effect 
on the PT03REA inner / PT03R outer combination. This 
result is surprising, considering that this combination has an 
entanglement consumption rate of 2 /4 while the entanglement 
consumption rate of the PT03R inner / PT03REA outer 
combination is smaller at 1/4. The result suggests that it would 
be wiser to place ebits in the inner encoder rather than in the 
outer encoder when these codes operate in practice. 

VII. Recursive, Classically-Enhanced Subsystem 
Encoders are Catastrophic 

We can construct other variations of quantum convolutonal 
encoders and study their state diagrams to determine their 
properties. One example of a variation mentioned at the end 
of Ref. |7 ] is a subsystem convolutional code (based on the 
idea of a subsystem quantum code |67j, (68), l69l that is 
useful in fault-tolerant quantum computation ITOll ). Poulin 
et al. suggest that subsystem convolutional codes might be 
"a concrete avenue" for circumventing the inability of a 
quantum convolutional encoder to be simultaneously recursive 
and non-catastrophic. Such codes exploit a resource known as 
a "gauge" qubit that can add extra degeneracy beyond that 
available in a standard stabilizer code. Another variation is 
an encoder that encodes both classical bits and qubits, and 
we might wonder if these could be simultaneously recursive 
and non-catastrophic (such codes are known as "classically- 
enhanced" codes im and are based on trade-off coding ideas 



from quantum Shannon theory [ 72 ]). 

Unfortunately, encoders that act on logical qubits, classical 
bits, ancilla qubits, and gauge qubits cannot possess both 
properties simultaneously, and we state this result below as 
a corollary of Theorem 1 in Ref. [7]. This result implies 
that entanglement is the resource enabling a convolutional 
encoder to be both recursive and non-catastrophic (there are 
no other known local resources for quantum codes besides 
ancilla qubits, classical bits, and gauge qubits). 

Corollary 1: Suppose that a classically-enhanced subsys- 
tem convolutional encoder is recursive. Then it is catastrophic. 

Proof: First, consider that the state diagram for a 
classically-enhanced subsystem quantum convolutional en- 
coder includes an edge from M to M' if there exists a 
/c g -qubit Pauli operator L q G {/, X, Y, Z} kq , an a-qubit 
Pauli operator S z G {/, Z} a , a k c -qubit Pauli operator 
L x + C z G {I,X,Y,Z} k \ and a #-qubit Pauli operator 
G x + G Z G {/, X, Y, Z} g such that 

(M' :P) = (M :L q :S z : L x + C z : G x + G Z )U, (5) 

where the binary representations of Pauli operators are the 
same as in ([3]). We break the operator acting on the k c 
classical bits into two parts because L x c represents the X 
component of the operator that can flip a classical bit from 
|0) to |1) and back, while C z represents the Z component 
of the operator that has no effect on a classical bit in state 
|0) or |1) (it merely adds an irrelevant global phase in the 
case that the bit is equal to |1)). The entries G x and G z are 
for the "gauge qubits" in a subsystem code [69] (qubits in 
a maximally mixed state that are invariant under the random 
application of an arbitrary Pauli operator). We include these 
transitions in the state diagram because a particular logical 
operator L = (L q : L x ) in a classically-enhanced subsystem 
code is equivalent up to Z operators acting on the ancilla 
qubits, Z operators acting on the classical bits, and X and Z 
operators acting on gauge qubits. The set C (L) corresponding 
to a particular logical input L = (L q : L x ) consists of all 
operators P of the following form: 

( P : P = (L q : S z : L x + C z : G x + G Z )V, \ 
C(L) = \ s z e{I,Z}\ C z e{I,Z} k % L 
[ G x + G Z e{I,X,Y,Z} 9 J 

where V in this case consists of the infinite, repeated, over- 
lapping application of the convolutional encoder U. 

Suppose that the encoder is recursive. By definition, it 
follows that every weight-one logical input L and its coset 
C (L) has an infinite response. Suppose now that we change 
the resources in the encoder so that all of the classical bits 
and gauge qubits become ancilla qubits. The resulting code is 
now a standard stabilizer code. Additionally, we can remove 
the X operators L x from the definition of C (L) because they 
no longer play a role as a logical operator, and we can remove 
the X operators G x from the definition of C (L) because they 
are acting on ancilla qubits. Let C' (V) denote the new set 
corresponding to a particular logical operator V = (L q : 

( P:P = (L q :S z :C Z : G Z )V, \ 
C'{L') = \ S z e{I,Z}\ C z e{I,Z} kc , \. 
{ G z e{I,Z} 9 J 
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The encoder is still recursive because C (L f ) C C (L f ) and 
because the original encoder is recursive (recursiveness is a 
property invariant under the replacement of cbits and gauge 
qubits with ancilla qubits). Then Theorem 1 of Ref. [7 ] implies 
that the encoder for the stabilizer code is catastrophic, i.e., 
it features a zero physical-weight cycle with non-zero logical 
weight. It immediately follows that the original encoder for the 
classically-enhanced subsystem code is catastrophic because 
its state diagram contains all of the edges of the stabilizer 
encoder, plus additional edges that correspond to the logical 
transitions for the classical bits and the X operators acting on 
the gauge qubits. ■ 

We should note that the above argument only holds if the 
original classically-enhanced subsystem convolutional encoder 
acts on a non-zero number of logical qubits (this of course is 
the case in which we are really interested in order to have 
a non-zero quantum communication rate). For example, the 
encoder in Figure [2] becomes recursive and non-catastrophic 
when replacing the logical qubit with a classical bit and the 
ebit with an ancilla. One can construct the state diagram from 
the specification in §5§ and discover that these two properties 
hold for the encoder acting on the cbit and ancilla. 

The above argument also does not apply if an ebit is 
available as an auxiliary resource for an encoder. We have 
found an example of a non-catastrophic, recursive encoder 
that acts on six memory qubits, one logical qubit, one ancilla 
qubit, one ebit, one cbit, and one gauge qubit. The seed 
transformation for this example is as follows: 

{2934387, 3034525, 2323986, 3804870, 964819, 2111270, 
3722470, 1017852, 2678995, 912985, 1575948, 3138834, 
142503, 317898, 3473283, 1492935, 2701400, 2492844, 
2078030, 1894781, 915544, 75029}. 

A non-catastrophic, non-recursive subsystem convolutional 
encoder with a high free distance could potentially serve well 
as an outer encoder of a quantum turbo code, but this might not 
necessarily be beneficial if our aim is to achieve the capacity 
of a quantum channel. The encoder for a subsystem code is 
effectively a noisy encoding map because it is a unitary acting 
on information qubits in a pure state, ancilla qubits in a pure 
state, and gauge qubits in the maximally mixed state. Ref. l52l 
proves that an isometric encoder is sufficient to attain capacity, 
and we can thus restrict our attention to subspace codes rather 
than subsystem codes. Though, this line of reasoning does not 
rule out the possibility that iterative decoding could somehow 
benefit from the extra degeneracy available in a subsystem 
code, but this extra degeneracy would increase the decoding 
time. 

VIII. Classically-Enhanced EAQ Encoders 

We might also wish to construct classically-enhanced EAQ 
turbo codes that transmit classical information in addition to 
quantum information, in an effort to reach the optimal trade-off 
rates from quantum Shannon theory 1521 , l53l . These codes are 
then based on the structure of the codes in Refs. ITU . 1731 . The 
state diagram for the encoder includes an edge from M to M' 
if there exists a k q -qubit Pauli operator L q G {/, X, Y, Z} kq , 
a /c c -qubit Pauli operator L x c G {/, X} kq , a /c c -qubit Pauli 



operator L z G {I,Z} kq , and an a-qubit Pauli operator S z G 
{/, Z} a such that 

(M f :P) = (M:L q :S z :L X C +L Z : I C )U. 

The logical operator L q acts on k q qubits while the operator 
L x + L z acts on k c classical bits. We break this latter operator 
into two parts for the same reason discussed in the previous 
section. The state diagram for an encoder of this form is 
similar to that for an EAQ convolutional encoder because it 
includes all transitions for the classical bits. The difference is 
in the interpretation of the logical weight of edge transitions 
and in the logical label of an edge. If an edge features a Z 
operator acting on a classical bit, then this operator does not 
contribute to the logical weight of the transition and does not 
appear on the logical labels — an / appears on the logical label 
if / or Z acts on the classical bit and an X appears on the 
logical label if an X or Y acts on it. 

We have found several examples of recursive, non- 
catastrophic encoders acting on these resources. One of our 
examples acts on five memory qubits, one logical qubit, 
one ancilla qubit, one ebit, and one classical bit. Its seed 
transformation is as follows: 

{162568, 174333, 260907, 261622, 168521, 3644, 
101236, 204504, 227699, 260074, 138781, 9252, 
257257, 179548, 78687, 6352, 230925, 180362}. 

IX. The Preservation of Recursiveness and 
Non-Catastrophicity under Resource 
Substitution 

The technique used to prove Corollary [T] motivates us to 
consider which resource substitutions preserve the properties 
of recursiveness and non-catastrophicity. We have two differ- 
ent cases: 

1) A resource substitution that removes edges from the state 
diagram preserves non-catastrophicity and recursiveness. 

2) A resource substitution that adds edges to the state 
diagram preserves catastrophicity and non-recursiveness. 

To see the first case for recursiveness preservation, consider 
a general encoder acting on logical qubits, ancilla qubits, ebits, 
cbits, and gauge qubits. Each logical input L = (L q : L x ) to 
the encoder has the following form: 

( P: 1 

rm= J P=(L q :S*:I c :LZ + C*:G* + G*)V, 

[ ) -\ S z e{I,Z}\ C z e{I,Z} k % [' 
[ G x + G z G {I,X,Y,Z} 9 

where the conventions are similar to what we had before. 
Suppose the encoder is recursive so that C (Xi : I), C (Y$ : I), 
C (Zi : I), and C (I : Xj) all have an infinite response. Then 
the encoder is still recursive if we replace an ancilla with 
an ebit because the original encoder is recursive and we no 
longer have to consider the Z operator acting on the replaced 
ancilla. The encoder is still recursive if we replace a cbit with 
an ancilla qubit because we no longer have to consider the 
coset C (I : Xj). Finally, it is still recursive if we replace a 
gauge qubit with an ancilla qubit because we no longer have 
to consider the X operator acting on the replaced gauge qubit. 
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To see the first case for non-catastrophicity preservation, 
suppose that an encoder is already non-catastrophic. Then 
a resource substitution that removes edges from the state 
diagram preserves non-catastrophicity because this removal 
cannot create a zero physical-weight cycle with non-zero 
logical weight. 

To see the second case for non-recursiveness preservation, 
suppose that an encoder acting on logical qubits and ebits 
is non-recursive, meaning that at least one of the weight- 
one logical inputs Xi, Yi, or Z{ has a finite response. Then 
replacing an ebit with an ancilla certainly cannot make the 
resulting encoder recursive because we have to consider the 
response to the operators Xi, Yi, or Zi and one of these is 
already finite from the assumption of non-recursiveness of 
the original encoder. Furthermore, consider a non-recursive 
encoder acting on logical qubits, ebits, and ancilla qubits. 
Replacing some of the ancilla qubits with ebits or gauge 
qubits cannot make the encoder become recursive for the same 
reasons. 

To see the second for catastrophicity preservation, suppose 
that an encoder is catastrophic. Then a resource substitution 
that adds edges to the state diagram preserves catastrophicity 
because any zero physical-weight cycles with non-zero logical 
weight are still part of the state diagram for the new encoder. 

The following diagram summarizes all of the above obser- 
vations. The resources are L for logical qubits, S for ancilla 
qubits, E for ebits, C for ebits, and G for gauge qubits. 
Recursiveness and non-catastrophicity preservation flow down- 
wards under the displayed resource substitutions, while non- 
recursiveness and catastrophicity preservation flow upwards 
(substitutions at the same level can go in any order). 



( 


L : 


: S 


: E : 


: C : 


: G ) 






It 




It 


It 


( 


L : 


: E 


: E : 


: S : 


: S ) 










It 


It 


( 


L : 


: E 


: E : 


: E : 


: E ) 



We can then understand the proof of Corollary [T] in the 
context of the above diagram. The original encoder acts on 
(L : S : C : G) and is assumed to be recursive. Resource 
substitution of C and G preserves recursiveness, while making 
the encoder act on only ancilla qubits (the auxiliary resource 
for a standard quantum code). Using the fact that standard 
recursive quantum encoders are catastrophic (7), we then back 
substitute the resources, which is a catastrophicity preserving 
substitution. 

X. Conclusion and Current Work 

We have constructed a theory of EAQ serial turbo coding 
as an extension of Poulin et a/.'s theory in Ref. (7). The 
introduction of shared entanglement simplifies the theory 
because an EAQ convolutional encoder can be both recursive 
and non-catastrophic. These two properties are essential for a 
quantum serial turbo code to have a minimum distance that 
grows near-linearly with the length of the code, while still 
performing well under iterative decoding. We provided many 
examples of EAQ convolutional encoders that satisfy both 
properties, and we detailed their parameters. We then showed 



how the concatenation of these encoders with some of Poulin 
et al.'s and some of our own lead to EAQ serial turbo codes 
with near-linear minimum distance scaling that outperforms 
the scaling of the codes in Ref. (74). We conducted sev- 
eral simulations of EAQ turbo codes — two notable surprises 
were that placing ebits in the inner encoder can achieve a 
better performance than expected in both scenarios with and 
without ebit noise. Our simulations are generally consistent 
with the findings in Ref. 1321 , [28] and other results from 
quantum Shannon theory (2T) . (50) . namely, that entanglement 
assistance can significantly enhance error correction ability. 
Finally, we considered how to construct the state diagram 
for encoders that derive from other existing extensions to 
the theory of quantum error correction, and we showed that 
classically-enhanced subsystem convolutional encoders cannot 
simultaneously be recursive and non-catastrophic. 

There are many questions to ask going forward from here. 
It would be interesting to explore the performance of the other 
suggested code structures in Section [VIlT| to determine if they 
could come close to achieving the optimal rates from quantum 
Shannon theory (52), (53) . For example, what is the best 
arrangement for a classically-enhanced EAQ code? Should 
we place the classical bits in the inner or outer encoder? 
Are there more clever ways to use entanglement so that we 
increase error-correcting ability while reducing entanglement 
consumption? Consider that Hsieh et al. recently constructed a 
class of entanglement-assisted codes that exploit one ebit and 
still have good performance (75) . 

We should stress that the behavior of the entanglement- 
assisted codes using maximal entanglement is exactly like 
that of a classical turbo code. There is no degeneracy, and 
the iterative decoding algorithm is exactly the same as the 
classical one. Furthermore, analyses of classical turbo codes 
should apply directly in these cases (4), (76) . and it would be 
good to determine the exact correspondence. These analyses 
studied the bit error rate rather than the word error rate, so any 
study of EAQ turbo codes would have to factor into account 
this difference. 

Much of the classical literature has focused on the choice 
of a practical interleaver rather than a random one (77), (78) , 
(79), (80) , and it might be interesting to import the knowledge 
discovered here to the quantum case. 

Finally, it would be great to find examples of EAQ turbo 
codes with a positive catalytic rate that outperform the turbo 
codes in Ref. (7), in the sense that they either have a higher 
catalytic rate while tolerating the same noise levels or they 
have the same catalytic rate while tolerating higher noise 
levels. We tried about five different combinations with no 
success for these examples. 
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